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The extinction efficiency for the interaction of a plane wave with a large nonabsorbing spherical particle is 
approximately 2.0. When a Gaussian beam of half-width wq is incident upon a spherical particle of radius 
a with wo/ a < 1, the extinction efficiency attains unexpectedly high or low values, contrary to intuitive 
expectations. The reason for this is associated with the so-called compensating term in the scattered field, 
which cancels the field of the Gaussian beam behind the particle, thereby producing the particle’s shadow. 
I introduce a decomposition of the total exterior field into incoming and outgoing portions that are free 
of compensating terms. It is then shown that a suitably defined interaction efficiency has the intuitively 
expected asymptotic values of 2.0 for Wo/a » 1 and 1.0 for wo/a « 1 . 


1. INTRODUCTION 

In Lorenz -Mie theory, extinction is the name given to the 
energy lost by a plane wave during its interaction with 
a single spherical particle. 1 * 2 The energy lost is said to 
be either absorbed or scattered. The absorbed energy ex- 
cites the internal structure of the particle. The scattered 
energy is carried away from the particle by the radially 
outgoing electromagnetic waves created by the interac- 
tion. The scattered waves are further subdivided into 
diffracted waves created by the interaction of the plane 
wave with the geometry of the particle (i.e., the shape of 
its projected area) and specularly reflected and transmit- 
ted waves created by the interaction with the particle’s 
composition (i.e., the strength of the interaction is deter- 
mined by the particle’s refractive index). If a plane wave 
is incident instead upon an ensemble of randomly posi- 
tioned particles, then, in the single-scattering approxi- 
mation, extinction describes the exponential attenuation 
of the undeflected portion of the plane wave as it passes 
through the ensemble. 3-5 

The purpose of this paper is to determine the physical 
interpretation of extinction if a Gaussian beam of half- 
width u>o is incident upon a single spherical particle of 
radius a. If the particle is nonabsorbing, we intuitively 
expect that the extinction efficiency (he., the extinc- 
tion divided by the incident energy striking the surface 
of the particle) should behave in the following way as the 
width of the beam is varied. For plane-wave incidence 
and in the large-particle limit a » A the extinction effi- 
ciency is approximately 6 * 7 2.0. Half of this value is due 
to deflection of the geometrical rays that strike the par- 
ticle’s surface, and half is due to diffraction of the rays 
that graze its edge. 9 For a narrow beam incident upon 
a large particle with w 0 « a the geometrical rays that 
strike the particle’s surface are again deflected. But 
since the portion of the beam that grazes the particle’s 
edge is exceedingly weak, diffraction is correspondingly 
weakened. Thus the extinction efficiency should ap- 
proach 1.0, being due solely to the deflected rays. When 
the extinction efficiency is computed for wq « a, this ex- 
pected reduction that is due to the weakening of diffrac- 
tion does not occur. Rather, the extinction efficiency as 


a function of the particle size parameter continues to os- 
cillate about 2.0. But the amplitude of the oscillations 
increases dramatically as wo/a decreases. 

The first goal of this paper is to help us to understand 
physically what extinction describes for Gaussian-beam 
scattering and why the amplitude of the oscillations in 
€« xt increases as wo/a decreases. The second goal is to 
construct another efficiency that has the intuitive limits 
of 2.0 for u>o » cl and 1.0 for u/q « a. Our pursuit of 
these goals proceeds in the following way. In Section 2 
we briefly review the derivation of the extinction efficiency 
for a focused Gaussian beam striking a spherical particle 
head on. We then compute the extinction efficiency and 
observe its behavior as a function of the particle size pa- 
rameter for various values of w Q /a. In Section 3 we show 
that this behavior results from associating extinction with 
the decomposition of the total electromagnetic fields ex- 
terior to the particle into a sum of two parts, the fo- 
cused Gaussian beam in the absence of the target particle 
plus the scattered field produced by the target’s presence. 
This decomposition is a natural choice when the scattered 
field only weakly perturbs the beam for wo/a » 1. But 
when wo/a < 1 and the particle blocks off a large por- 
tion of the beam, the scattered field is no longer a weak 
perturbation, and this decomposition leads to certain diffi- 
culties in interpretation. In Section 4 we make an alter- 
native decomposition of the electromagnetic fields into a 
sum of two other parts, a radially incoming part and a ra- 
dially outgoing part. We define the interaction efficiency 
to be the energy carried by that portion of the outgoing 
wave that is distinguishable from the undeflected portion 
of the beam, divided by the incoming energy that strikes 
the particle’s surface. We show that the interaction effi- 
ciency behaves in a way consistent with our intuition in 
both the Wo/a » 1 and the w 0 /a « 1 limits. 

2. EXTINCTION EFFICIENCY FOR 
GAUSSIAN-BEAM SCATTERING 

A. Beam When No Target Particle Is Present 
Consider a focused Gaussian beam propagating along the 
z axis of a coordinate system. This is called an on-axis 
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beam. The time dependence of its electric and magnetic 
fields is exp(—iojf)i which will hereafter be omitted. The 
beam has wavelength A. It is focused by a lens to the 
half-width wq at the origin. One description of the elec- 
tric and magnetic fields of the beam is given by the Davis 
first-order approximation 9 


So exp (ikz) 
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The plane-wave limit is obtained by setting w o — * x and 
s = 0 in Eqs. (1). 

The beam fields have the partial- wave expansion 10 
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In these expressions the ji(kr) are spherical Bessel func- 
tions, and the angular functions irt(9) and rdd) are related 
to associated Legendre polynomials by 


ir t (0) 


Pi l { cos 6) 
sin 9 
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The shape of the beam is determined by the coefficients 
gi, which are weighting factors for the individual partial 


waves. Alternatively, these beam shape coefficients may 
be obtained from the beam profiles by 11 
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x exp(iAr cos ^)P/ l (cos 9), (6) 

where the radial components of the electric and magnetic 
fields of the on-axis beam assume the forms 
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For an on-axis Gaussian beam focused at the origin the 
beam shape coefficients are approximated to a high degree 
of accuracy 11 by the localized beam model 12 

gi * exp[-s 2 (/ + 1/2) 2 ]. (8) 

The partial-wave expansion of a plane wave is obtained 
in the limit s = 0 and gi = 1. 

B. Scattered Wave Produced by the 
Beam-Particle Interaction 

Consider now that the Gaussian beam approximated by 
Eqs. (1) or relation (8) substituted into Eqs. (4) is inci- 
dent upon a spherical particle of radius a and refrac- 
tive index n whose center is at the origin of coordinates. 
The standard method for solving for this electromagnetic 
boundary-value problem begins with the decomposition of 
the electric and magnetic fields exterior to the particle 
into a sum of two parts 1 *" 15 : 


E total — Ebeam + E Ka ttered » Btotal B^eam Blattered • 

(9) 


The fields Ebeem and Bbeam are the electromagnetic fields 
of the beam in the absence of the particle as given in 
Eqs. (4), and E^neied and Bettered represent the change 
in the total fields that is due to placing the target particle 
in the beam’s path. Specifically, the decomposition of 
Eqs. (9) gives 10 
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where the h\ l \kr) are spherical Hankel functions of the 
first type and a t and b[ are the partial-wave scattering 
amplitudes of plane-wave Lorenz- Mie theory. 

The total exterior field is illustrated in Fig. 1 for the 
case in which the beam is wide and the particle is small. 
At most locations in space the total fields strongly re- 
semble the original beam, and B^ ttere< i repre- 

sent only a small perturbation, and the decomposition of 
Eqs. (9) makes good physical sense. An exception, how- 
ever, is provided by the deep-shadow region immediately 
behind the particle, where the particle blocks off the ori- 
ginal beam. 16,17 The deep-shadow region is denoted by S 
in Fig. 1. In spite of the particle's blocking off the inci- 
dent beam in the deep-shadow region, the total field there 
is not necessarily small if the particle is transparent. In 
particular, surface waves associated with reflection pass 
through this region (see Fig. 46 of Ref. 8), and the trans- 
mitted waves form a spherical aberration caustic there. 
The term deep-shadow region refers instead to the fact 
that the waves diffracted by the particle only gradually 
spread to fill the space behind it (see Fig. 14 of Ref. 17). 
Similarly, the term block off means that, although the 
fields in the deep-shadow region are nonzero, they do not 
at all resemble the incident beam. The incident beam 
has been totally removed and replaced by the reflection 
surface waves and transmitted waves. As a result, in the 
deep-shadow region the scattered fields must assume the 
form 16 


E scattered ** E com p ensa ting + s ,d U »l 
Bscattered ~ Bcompensating + Bresidual 


(ID 


where 


Ecomp«niating ” Ebetm » Bcompcnaating * -b„, am • ( 12 ) 

The compensating term in the scattered field cancels the 
beam field and thereby mathematically ensures its re- 
moval from the deep-shadow region. The residual fields 
are the transmitted waves and the reflection surface 
waves present in the deep-shadow region if the particle is 
transparent. If it is opaque, the residual fields are only 
the reflection surface waves. 17 


C. Extinction Efficiency 

Using the decomposition of Eqs. (9), we can define the ab- 
sorption, scattering, and extinction cross sections as 2 * 10 18 
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respectively. They reduce to the usual Lorenz -Mie ex- 
pressions in the plane-wave limit gi = 1. Since the 
beam's angle-integrated Poynting vector vanishes, i.e., 

f W sin Odd (*' d<t> Re(E^ ara x = 0. (16) 

Jo Jo 

the absorption, scattering, and extinction cross sections 
are related to one another by 

^extinction “ (-^scattering (^absorption * (1/ 

Equation (17) describes conservation of energy for the 
beam -particle system. The absorption, scattering, 
and extinction efficiencies are obtained by dividing 



Fig. 1. Scattering of a plane wave by a spherical particle 
The small region S to the right of the particle denotes the 
deep-shadow region, where the plane wave has been removed 
by the particle in the near-forward direction and in the near 
zone and replaced by the transmitted wave and the reflection 
surface waves. 
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Eqs. (13)— (15) by the cross section for radiation strik- 
ing the particle surface, 
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In the Davis first-order beam approximation the incident 
cross section is 
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In the plane- wave limit s — » 0, Eq. (20) reduces to 

^incident * 

The behavior of the extinction efficiency of Eq. (19) as 
the width of the beam is varied is illustrated in the fol- 
lowing example. Consider a spherical water droplet of 
radius a * 50 fixti and n ** 1.333 illuminated on axis by a 
Gaussian beam focused at the origin with 0.57 /xm ^ A ^ 
0.64 jim and 0.1 ^ w Q /a ^ 100. These values of w 0 are 
well within the region of applicability of both the Davis 
first-order beam approximation of Eqs. (1) and the local- 
ized model for the beam shape coefficients of relation (8). 


The incident cross section of Eq. (20) was computed by nu- 
merical integration with a 5001 -point grid. The result- 
ing extinction efficiency of Eq. (19) is shown as a function 
of the particle size parameter 


for wo/a « 20.0, 1.0, 0.2, and 0.1 in Figs. 2(a), 2(b), 2(c), 
and 2(d), respectively, and as a function of wo/a for 
A - 0.6328 /xm in Fig. 3. For w Q /a » 1 we find that 
6„ t =*• 2.0, in agreement with the plane-wave limit. This 
result is interpreted as being an amount 1.0 that is due 
to scattering plus another 1.0 that is due to diffraction. 6 
In addition, the extinction efficiency as a function of size 
parameter possesses oscillations, known as the inter- 
ference structure, that are due to interference between 
the diffracted and transmitted fields. 19,20 The efficiency 
also exhibits the so-called ripple structure that is due to 
morphology-dependent resonances. 21 For w Q /a < 0.5 the 
extinction efficiency as a function of size parameter con- 
tinues to oscillate about 2.0. But the amplitude of the in- 
terference structure increases dramatically. The ripple 
structure is now absent, since the Gaussian beam no 
longer extends out to the edges of the particle, where 
morphology-dependent resonances are most efficiently 
excited. 22 The increase in the amplitude of the inter- 
ference structure is contrary to our expectation that the 
diffraction contribution to the extinction efficiency should 
decrease because a progressively weaker portion of the 
beam grazes the edge of the particle. In Section 3 below 
we demonstrate that the large-amplitude interference 
structure for wo/a « 1 results from the interference of 
the transmitted field with the compensating field. 
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W 0 /■ 

Fig. 3. Extinction efficiency of Eq. (19) as a function of the 
half-width of the beam divided by the particle radius for a beam 
wavelength of A * 0.6328 /im, a particle radius of a — 50 /im, 
and a refractive index of n - 1.333. 



Fig. 4. Scattering of a narrow beam by a spherical particle. 
The region B to the right of the particle and extending to the far 
zone denotes the region where the beam has been removed by the 
particle and replaced by the transmitted and reflected waves. 


3. COMPENSATING FIELD IN 
GAUSSIAN-BEAM SCATTERING 

Figure 4 illustrates the interaction of a narrow beam with 
a large particle. The particle blocks off most of the beam, 
preventing its continued propagation to z — - «. Again the 
term block off means that the beam behind the particle 
has been removed and replaced by the transmitted and 
reflected waves. On the other hand, the decomposition 
of the total fields in Eqs. (9) contains Eb*an, and Bbe*in, 
which in fact are not present behind the particle. Thus 
E scattered and Blattered in the near- forward direction and in 
both the near zone and the far zone must be of the form 
of Eqs. (11) so that the compensating field can cancel the 
beam field and thereby ensure the beam's absence behind 
the particle. 16 

The compensating field is easily identified in numerical 
computations of the far-zone scattered intensity. Con- 
sider a focused on- axis Gaussian beam with u> 0 = 10 /xm 
and A = 0.6204 /im incident upon a spherical water drop- 
let with a «* 50 /xm and n = 1.333. This wavelength is 
equivalent to the size parameter x =* 506.38 and corre- 
sponds to the first maximum of the extinction efficiency 
graph of Fig. 2(c), where - 2.734. The far-zone scat- 


tered intensity 

Ixmtt ered(0> = lim X Bettered) (22) 


is graphed in Fig. 5(a) for - 180* 5 S < 180* and <t> = 
90*. The scattering is dominated in the forward hemi- 
sphere by transmission and in the backward hemisphere 
by specular reflection and transmission following one in- 
ternal reflection. 11 Since w 0 /a « 1, diffraction in the 
near-forward direction is minimal. In Fig. 5(a) there is 
a large peak at -1.5* ^ 9 1.5* that is not observed in 

experiments. 23 * 24 I claim that this peak is the intensity 
corresponding to the compensating field. The inset of 
Fig. 5(a) shows the scattered intensity for 0* < 6 < 5*, il- 
lustrating the constructive interference between the com- 
pensating field (C) and the transmitted field (T) for A = 
0.6204 /im. Since the scattered intensity for -5* < 6 < 
5* is orders of magnitude larger than at any other scatter- 
ing angle, integrating over this constructive interference 
produces the increased extinction efficiency = 2.734 
shown in Fig. 2(c). 

Similarly, the far-zone scattered intensity for A = 
0.6261 /xm is shown in Fig. 5(b). This wavelength is 
equivalent to x — 501.77 and corresponds to the first 
minimum of the extinction efficiency graph of Fig. 2(c), 
where e txt = 1-271. The inset of Fig. 5(b) illustrates 
the destructive interference between the compensating 
field (C) and the transmitted field (T) for A = 0.6261 /xm. 
Integrating over this destructive interference produces 
the decreased extinction efficiency of e <xt = 1-271. A de- 
tailed physical model of the oscillations in the extinction 
efficiency as a function of size parameter for w Q /a « 1 
based on compensating field- transmission field interfer- 
ence will be given elsewhere. 

As interesting as these large-amplitude oscillations in 
for wq/q « 1 appear to be, they are not observable. 
Figures 5(c) and 5(d) graph the far-zone intensity associ- 
ated with the total field exterior to the particle, 25 

= lim RefE^i x B total ), (23) 


for -180* ^ S ^ 180* and <t> = 90* for A =* 0.6204 /xm and 
A = 0.6261 /xm, respectively. Figures 5(c) and 5(d) are 
identical to Figs. 5(a) and 5(b) except that the beam fields 
have canceled the compensating portion of the scattered 
fields and have eliminated the spurious peak in Figs. 5(a) 
and 5(b). This cancellation verifies my earlier claim of 
having properly identified the compensating field peak in 
Figs. 5(a) and 5(b). It also affects the physical meaning 
of the extinction efficiency. Since the large-amplitude in- 
terference structure in e ext was produced by the interfer- 
ence between the compensating field and the transmitted 
field and since only the total field is observed in experi- 
ments, the compensating field- transmission field inter- 
ference is not observable. If another efficiency based on 
the total field were to be constructed, the cancellation of 
the beam by the compensating field would prevent the 
large-amplitude interference structure in the other effi- 
ciency from occurring. 

The compensating term in the scattered field in the 
near-forward direction for wo/a s 1 can also be demon- 
strated analytically. For small S the angular functions 
of Eqs. (5) are approximated by 26 


*r,(0) * r,(0) 
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expansion of the plane- wave Mie theory partial- wave scat- 
tering amplitudes is 28 * 29 


ai 

bi 
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where the first term (1) denotes diffraction, the second 
term (-Rf 2 ) denotes specular reflection, and the third 
term [-Tf 1 (i2/ 1 ) P_1 T/ 2 ] denotes transmission following 
p-l internal reflections. The diffracted portion of the 
scattered field is then 

^Scattered * ^ exp(ikr)u x £ (Z + 1/2 )gtJo[(l + 1/2)0]. 

* r Z-l 

(34) 


Again converting the sum over partial waves to an inte- 
gral over an associated impact parameter u, we obtain 


general present at a given point in space, they are distin- 
guishable because the wave fronts of the outgoing beam 
are flat and its amplitude is constant, whereas the wave 
fronts of the scattered wave are spherical and its ampli- 
tude falls as 1/r (see Fig. 1). For a narrow beam incident 
upon a large particle with wja « 1 the beam and the 
scattered parts are not individually observable, since the 
particle has removed the beam in the near-forward di- 
rection and replaced it by the transmitted and reflected 
waves. The removal of the beam is not evidenced in ei- 
ther El*.*,, or Entered taken individually but is seen only 
when they are added together to form E^i. In experi- 
ments, however, only the total field is observed. The 
extinction efficiency exhibits interference between the 
compensating and transmitted portions of the scat- 
tered field that cannot be observed because of the beam 
field -compensating field cancellation. 
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diff 
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where 30-32 
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Since the spherical Bessel function ji(kr) in Eqs. (4) may 
be written as 
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Jo 


wq > a 


1 ,2/ V f 2 a 2 V J n (kad) 

— exp(-aW)I^J 
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For a wide beam incident upon a small particle with 
wq/cl » 1 the diffracted electric field is approximately 

Eti? tter «i * exp(i kr)u x exp(-a 2 /wo 2 )(ka) 2 ' 

(37) 

This is the usual state of affairs in plane-wave Mie 
theory. But for wq « a a comparison of relation (28) 
with relation (35) and Eq. (36) leads us to identify the 
first term in the second line of Eq. (36) substituted into 
relation (35) as the compensating field. After the can- 
cellation of the compensating field and the beam field the 
remaining part of Ei£ Ured , given by 

E*c»tured md * r “ exp (ikr)u x exp(-a 2 /u> 0 2 ) 

x J 0 (ka6) , (38) 

describes the weak diffraction produced by the tail of the 
beam grazing the edge of the particle. A qualitative and 
intuitive derivation of relation (38) is given in Appendix A 
below. 

It is now clear what extinction in Gaussian-beam scat- 
tering describes and why the extinction efficiency of 
Eq. (19) behaves in a counterintuitive way for wja < 1. 
Extinction is based on the decomposition of the total fields 
exterior to the particle into the sum of a beam part and a 
scattered part as in Eqs. (9) and on the implicit assump- v 
tion that each part individually is observable. For a 
wide beam incident upon a small particle with wja » 1 
this is indeed the case. Although both the undeflected 
portion of the original beam and the scattered wave are in 


ji(kr) = l f 2 h\ l \kr) -r l / 2 hf\kr ) , (39) 

where h[ l) (kr) describes radially outgoing waves and 
h! 2 \kr) describes radially incoming waves, 33 by substi- 
tuting Eq. (39) into Eqs. (4) we may decompose the total 
field into the sum of an incoming part and an outgoing 
part, 

^ total = ' ^outgoing + E mcomin g , 

®totai ^outgoing +* Bincoming • (40) 

The fields E,*,^,,^ and B outgoing contain all the h^ikr) 
terms, and the fields E mcomin g and Burning contain all 
the h i2 \kr) terms. The partial-wave expansion of the 
outgoing fields is identical to Eqs. (10) except that a* and 
bi are replaced by 

A/* a t -1/2, B, 1/2, (41) 

respectively. Similarly, the partial-wave expansion of 
the incoming fields is identical to Eqs. (4) except that 
ji(kr) is replaced by l hh\ 2) (kr). 

In a number of respects the incoming and outgoing 
fields have a pleasing physical interpretation. First, the 
Debye-series expansion of the incoming and outgoing 
partial-wave amplitudes is sensible. For the incoming 
fields, using a shorthand notation, we have 

e ~— <421 

In the large-particle limit relation (42) reflects the fact 
that the incoming beam consists of geometrical light rays 
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that will strike the particle (first term) and rays that will the incoming, outgoing, and absorption cross sections are 

miss it (second term). For the outgoing fields we have related by 


E, 


outgoing 




Relation (43) reflects the fact that the outgoing radiation 
consists of geometrical rays that have struck the particle 
and have been either reflected or transmitted following 
p - 1 internal reflections (first term), rays that were 
incident at the edge of the particle and participated in 
tunneling reflection or tunneling transmission following 
p - 1 internal reflections 3 * 34 (second term), and diffracted 
rays plus the outgoing portion of the original beam that 
missed the particle (third term). In the large- particle 
limit the third term in relation (43) may be written as 

E^oSr 1 = T 1 exp (ikr)u z £ (l + 1/2)*, (-1) 

Kr l-tmi 

x J 0 [(l + 1/2)0] 

= - exp(i£r)u x f udu exp(-s 2 u 2 )Jo(u0) 

kr Jka 

=** expfifcr)^ j udu exp(-s 2 u 2 ) 

f ka <> 

X J o (u0) - / udu exp (-s 2 u 2 )J o (u0) . 


For wo/a » 1 the first term of the last two lines of 
Eq. (44) is the undeflected portion of the beam [see rela- 
tion (28)], and the second term is the circular-aperture 
Fraunhofer diffraction field of relation (37). For wa/a < 
< 1 the last two lines of Eq. (44) contains only the remain- 
der diffraction term of relation (38). The outgoing fields 
in relation (43) are thus free of the unobservable compen- 
sating term that was a necessary part of the scattered 
fields for wja < 1 in the decomposition of Eqs. (9). 

The decomposition of the total fields into a sum of 
incoming and outgoing waves also provides a sensible 
description of conservation of energy for scattering by a 
single particle. Let us define the incoming and outgoing 
cross sections as 

c — « = 0/ o sm dds L H 

x RaCE*^ X B incoming) 

= TT i(i + V2)l*/l 2 , (45) 

* /-I 

C ou ^ f --^/ o ' S in edef* d* 
x Re(E 0utf0ing x Boutgoing) 

-T7 i(2/ + l)|*,lW + l5,l 2 ), (46) 

* 2 it] 

respectively. Since the angle-integrated cross-term 
Poynting vector vanishes, i.e., 

ttt f sin ddd f d <t> 
to 1 Jo Jo 

X Re (E ^coming X Boutgoing + E outgomg X B incoming) = 0, (47) 


c 


incoming 


c 


outgoing 


+ c 


absorption ■ 


(48) 


The outgoing efficiency is defined as 


c outgoing 


outgoing 

incident 


c, 


(49) 


where is given by Eqs. (18) and (20). 

The outgoing efficiency of Eq. (49) has a problemati- 
cal behavior for wo » a. But a redefinition of ^outgoing 
along the lines suggested more than 40 years ago by 
Brillouin 16 removes the problem. The problematical be- 
havior of ^outgoing may be easily demonstrated by the 
following example. Equation (49) was calculated for a fo- 
cused on-axis Gaussian beam with A = 0.6328 pm and 
0.1 ^ wo/ a ^ 10 incident upon a spherical water droplet 
with a m 50 pm and n =* 1.333. The results are shown 
in Fig. 6. For wo/a < 1 the outgoing efficiency is 1.0, 
being due almost entirely to geometrical rays striking the 
particle surface and deflected by it. Diffraction is negli- 
gible, For wo/ a ** 1 the beam at the edge of the particle 
is stronger, diffraction becomes correspondingly stronger, 
and e 0 utgoujg increases. This is sensible behavior. But 
for wo/a > 1 the outgoing efficiency increases beyond 
2.0, in disagreement with one's intuitive expectation of a 
maximum value of approximately 2.0. 

The problem is that for wo/a > 1, in the near-forward 
direction, the outgoing fields of relation (43) contain not 
only the diffracted, reflected, and transmitted fields but 
also the portion of the original beam that missed the par- 
ticle and now is propagating undeflected toward z * as 
in Fig. 1. If we wish to describe the efficiency of the in- 
teraction of the beam with the particle, this undeflected 
field must be subtracted from E 0utg0l ng. I propose to do 
this in the following way. The far-zone outgoing inten- 
sity is constructed from the outgoing fields according to 


foutgoing(0r “ ^(ISoutgomg X B ou tgoing) • (50) 


Since the two largest contributions to the outgoing in- 
tensity in the near-forward direction for wo/a » 1 are 



Fig. 6. Outgoing efficiency of Eq. (49) as a function of the 
half-width of the beam divided by the particle radius. The beam 
and particle parameters are the same as those in Fig. 3. 
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(a) 



Fig. 7. (a) Angle of the first minimum of the outgoing intensity 

of Eq. (50) as a function of the half- width of the beam divided 
by the particle radius (solid curve). The beam and particle 
parameters are the same as those in Fig. 3. The dashed curve 
is the approximation of relation (53). (b) Interaction efficiency 

of Eq. (52) as a function of the half- width of the beam divided by 
the particle radius. The beam and particle parameters are the 
same as those in Fig. 3. 


the undeflected original beam of relation (28) and the dif- 
fracted wave of relation (37), the first relative minimum 
of /outgoing as a function of 6 y which we term 0 min , results 
from the destructive interference of the undeflected beam 
with the diffracted field and is nearly independent of 4>. 
Following Ref. 16, we claim that any scattering contri- 
bution to the outgoing intensity in the angular region 
0 * ^ 0 ^ 0 min and 0 s <f> ^ 2n cannot be disentangled 
from the undeflected original beam, since they both de- 
crease as 1/r. On the other hand, the outgoing intensity 
for 0 m in ^ 6 ^ 180* is easily recognized as being due to 
scattering, since it lies beyond the angular interval within 
which the undeflected beam is confined. When the out- 
going intensity is integrated over this range of 6 and <f > , 
the result represents the energy lost from the forward- 
propagating and spreading Gaussian beam during its in- 
teraction with the particle. We define this interaction 
cross section and efficiency to be 


Gaussian beam with A ==* 0.6328 /xm and 0.1 < wo/a < 
100 incident upon a spherical water droplet with a =* 
50 fji m and n * 1.333. The angle 0 min was numerically 
determined and is shown as a function of wo/a in Fig. 7(a). 
For w Q /a » 1, 0 min may be approximately obtained by 
locating the angle of maximal destructive interference of 
the undeflected original beam of relation (28) with the 
diffracted field of relation (37). The result, 



exp(-a 2 /wo 2 ) 


(53) 


is shown as the dashed curve in Fig. 7(a). As diffraction 
weakens for wq/cl < 1, the value of 8 min is determined 
by the interference between the original beam and the 
transmitted field, which is the next-largest contribution to 
E outgoing* For wo/a « 1, 0 min decreases to zero at the size 
parameters for which the beam field and the transmitted 
field destructively interfere as in Fig. 7(a). On the other 
hand, at the size parameters for which the beam field and 
the transmitted field constructively interfere for wo/a « 
1, the value of 0 mm levels out at a small nonzero value. 

The interaction efficiency of Eq. (52) was computed 
with a 7200-point grid for the 8 integration in extraction- 
The results are shown in Fig. 7(b). For w 0 /a » 1 
the interaction efficiency is approximately 2.0. Of this 
amount, approximately 1.0 is due to scattering and ap- 
proximately 1.0 is due to diffraction. For w 0 /a *= 1 
diffraction weakens, and the interaction efficiency corre- 
spondingly decreases. For wo/a « 1 virtually no diffrac- 
tion occurs, and the interaction efficiency is due almost 
entirely to scattering. Similar results occur for other 
wavelengths and particle sizes and are consistent with 
our intuition about diffraction for both wide and narrow 
beams. I claim that, for Gaussian-beam scattering, the 
interaction cross section of Eq. (51) represents the best 
measure of the energy lost by the incident beam resulting 
from its interaction with a single spherical particle. 

This entire development was for an on-axis Gaussian 
beam. The extinction cross section for an off-axis Gauss- 
ian beam has been derived in Ref. 35. Although we have 
not numerically analyzed this case, the extinction effi- 
ciency presumably also behaves in a counterintuitive way 
for Wo/a < 1, and a corresponding interaction efficiency 
may also be defined. 


APPENDIX A 

Consider diffraction of an incident beam by an aperture 
in the Fresnel -Kirchhoff approximation 36 

E diff (z) - ^ j dVEi» e jd.m(r')A(r')exp(-i^Ur • r'), 

(Al) 


C 


interaction 


^interaction 


_ * f sin Odd f d<£ /outing (0i <£)» 

&0 Jo 

(51) 


'-'interaction 

^incident 


(52) 


This subtraction procedure for the undeflected field was 
numerically tested in the following way. The outgoing 
intensity of Eq. (50) was computed for a focused on-axis 


where r is the vector from the center of the aperture plane 
to the position of the observer a distance z away, r' is the 
vector from the center of the aperture plane to any point 
in the aperture, E inC j dent (r / ) is the incident electric field 
in the aperture plane, and the aperture function A(r') is 
defined by 


A(r') 


1 inside the aperture 
0 outside the aperture 


(A2) 
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Let the incident beam be a focused on- axis Gaussian beam 
whose electric field in the aperture plane is 

Emcid.nt(r) - Eo exp(-r' 2 / wq 2 )u x exp (ikz), (A3) 

and let the aperture be the region exterior to a circle of 
radius a. This is appropriate for describing diffraction by 
a spherical obstacle. If the observer is at the scattering 
angles (0, 4>) t Eq. (Al) becomes 

E diff (z) - f’r'dr' exp(-r' 2 /wo 2 ) 

2ttz J o Jo 

x exp[-ifcr' sin 0 cos(f - exp(ifcz) 

=* u x f r'dr' exp(-r /2 /w ; o 2 Vo(^ / sin 0) 

Z Ja 

x exp (ikz). (A4) 


The strength of the incident electric field at the edge of the 
aperture is exp(-a 2 /u>o 2 ). It falls to 1/e of this value at 

r' * a + Ar' *■ a 4- ~~~ * (A5) 

The value of the last integral of Eq. (A4) is then crudely 
approximated by 


J r’dr 1 exp(-r' 2 /u>o 2 )«/o(^ sin 0) 


“( 


u v 
2a 


\exp(-a 2 /wo 2 )Jo(kad ) . (A6) 


giving 


E diff (z) 


exp(t&z) ^ exp(-a 2 /u/o 2 )*/o(^a0)& x » 

kz 2 

(A7) 


in agreement with relation (38). 
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